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This is a summary of the latest results of the Delphi Collaboration on the properties of identified 
quark and gluon jets. It covers the measurement of the fragmentation functions of giuons and 
quarks and their scaling violation behaviour as well as an analysis of the scale dependence of the 
multiplicities in gluon and quark jets. Further, a precision measurement of Ca/Cf from the 
multiplicities in symmetric three jet events is discussed. 



I. INTRODUCTION 

o\ : 

In QCD, the three fundamental splittings of the two types of colour charged fields participating in the strong 
interaction (quarks (q) and giuons (g)) are q — > qg, g — > gg, and g — > qq. The corresponding splitting kernels, which 
describe both the kinematics and the relative strengths of these splittings, are proportional to the colour factors 
Cf = § , C a = 3, and Tr = \uf, respectively, where Up is the number of active quark flavours in g — > qq decays. The 
values for these colour factors originate directly from the group SU(3) underlying QCD so that precise measurements of 
these couplings validate QCD as the fundamental theory for strong interactions. In the sensible range of uf = 3 ... 5, 
the ratio of the splitting kernels S( g ^, gg+g ^ qq -)/ S( q ^ qg ) is nearly constant (2.0 - 2.3) and exactly Ca/Cf =2.25 in limit 
j of large energy transfers. This leads to the expectation of a Ca/Cf times higher probability of gluon bremsstrahlung 
0^ ' in gluon jets with respect to quark jets. This ratio should be visible for observables which are proportional to the 
splitting probabilities of giuons and quarks 

Unfortunately, quarks and giuons are no free particles. Therefore one has to use jets of giuons and quarks in three 
jet events as the best approximation of the tree level graphs. Beyond the necessity to rely on the assumption of LPHD 
(Local Hadron Parton Duality) which states that the parton properties are preserved to the hadron level which is 
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then clustered to jets, this approach is further limited by the influence of interference effects in the event. Moreover 
the assignment of particles to jets is somewhat arbitrarily. 



' Further, the jet finding algorithms introduce ambiguities when assigning particles to jets 

43 : 



II. EXPERIMENTAL ACCESS TO GLUON AND QUARK JETS 



■ Gluon jets were originally identified in symmetric three jet events ( Y 




events) |2h4]. In these events the most energetic jet is excluded from the 
3 analysis, as the rate of gluon induced leading jets is rather low. Due to the 
symmetric event topology (#2 ~ #3, see Fig. 1 ), the low energy jets are 
:t2 expected to be directly comparable. By identifying one of these jets as a b- 

jei L 

quark jet using impact parameter techniques, the remaining jet is identified 

-^ T „ , Yeve "i s qqg events indirectly as a gluon jet. The properties of a comparable ucisc-quark jet 

FIG. 1. Symmetric and asymmetric J b J 

event topologies sample can then be obtained from non-6 events, where the gluon properties 

are eliminated from by subtraction techniques 
Recently this technique has also been extended to non-symmetric events Q] (see Fig. 1 ). In this case one relies on 
the quark/gluon composition of jets as predicted by the three jet matrix elements taken from Monte Carlo simulations. 
This technique improves the available statistic and gives access to a wider range of energy scales, but requires a criterion 
for the selection of comparable gluon and quark jets. These jets are obtained from events with different topologies. 
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The comparison of jet properties obtained for this event sample^ and for Y events helps in finding a suitable scale 
to classify the jets with. One yields about 20,000 identified gluons in Y events and about 100,000 in the asymmetric 
events. 

III. JET SCALES 

The relevant scale for the jet evolution is not just the jet energy. In Fig. ||, the fragmentation functions of quark jets 
as a function of their energy are shown. The rows in this figure correspond to data taken in the same xe intervals. The 
much more pronounced scaling violation pattern of jets obtained from symmetric events (Ej et ~ 22GeV . . . 29GeV) 
compared to jets from asymmetric events clearly disfavours the jet energy as the relevant scale. 
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FIG. 2. Quark frag. func. vs. Ej e 
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FIG. 3. Quark frag. func. vs. Kg 
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Soft radiation is limited to cones given by the opening angles between the jets. This motivates transverse momentum 
scales. The hardness kh = Esm6i/2 is a better choice ||, as it accounts for the limited phase space available for 
gluon radiation due to the interference of radiated gluons in the eventQ. Fig. || shows the good agreement of light 
quark jets in three jet events with jets in e + e~ annihilations from Petra energies || to the highest LEP energies (Tc| l 
in normalization and slopes using kh for the first and E\, eam for the latter. This agreement confirms the interpretation 
of kh as a valid scale for jet evolution in three jet events. In multi-jet events several scales may be relevant; in so 
far the usage of kh as a (single) scale is an approximation. Another possible choice is the scale pf , introduced 



in Sec. IV E 



IV. RESULTS 



Recently, several properties of quark and gluon jets which can be expected to be sensitive to Ca/Cf , such as the 
scaling violation of quark and gluon jets |Tl| and the scale dependence of the quark and gluon jet multiplicities E^] 
have been analysed using data collected with the Eelphi detector at LEP. 



1 In the following called asymmetric event sample though containing all Y events as well. 

This definition corresponds to the beam energy of an e + e~ — > qq event. Often also the definition k = 2i5sin(9/2 ~ E9 is 
used. 
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A. Gluon and Quark Fragmentation Functions 



The gluon and quark fragmentation functions D g>q (xE, kh) are measured in dependence of different values of kh 



at a fixed CMS energy (see Sec. Ill for a short discussion of jet scales) 



Y topology 
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FIG. 4. Gluon fragmentation functions for different values 
of Kg 
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FIG. 5. Gluon and quark fragmentation function at a fixed 
topology, 6*2 , 83 G [150° ± 15°] ), compared to different Monte 
CArlo generators 



The predictions are derived from the numerical solution of the DGLAP evolution equation in first order || . The 
following ansatz has been used to parameterize the fragmentation function at a fixed scale kq to start the evolution: 



D g p ' q (x E ) = P 9 3 ' q ■ xf ■ (1 - x E fi q ■ exp (~pf q ■ In 2 x E ) . 



(1) 



The parameters pf g ', Aqcd and the colour factor Ca are fitted simultaneously. 

The softening of the fragmentation functions with increasing ft is observed. This effect is more pronounced for 
gluon jets than for quark jets. Fig. || compares the gluon and the quark fragmentation function at a fixed topology. 



B. Scaling violation of the gluon and quark fragmentation functions 

Fig. H shows the measurement for light quark jets in three jet events as a function of kh, represented by the open 
symbols and superimposed by power law fits. The corresponding result for gluons is shown in Fig. ^ for symmetric 
and asymmetric event topologies. The good agreement between the two event samples in Fig. ^ demonstrates that 
Kh is a sensible choice for the scale of gluon jets as well. 

A simultaneous fit with the first order DGLAP equations to both the quark and gluon fragmentation functions has 
been performed. Beyond the parameters of the analytic ansatz of the fragmentation functions , Aqcd and Ca have 
been treated as free parameters. The fit is sensitive to the occurrence of Ca in the g — > gg splitting kernel. The fit 
describes the data well and yields: 
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FIG. 6. Gluon fragmentation functions as a function of jet 
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FIG. 7. (Logarithmic) slopes of the scale (kh) dependence 



scale kh ■ Superimposed are jets from both symmetric and of the quark and gluon fragmentation functions as a function 
asymmetric events of xe 

= 2.44 ± 0.21 sfat (preliminary) . 

Cf 

for the colour factor ratio. 

In Fig. the logarithmic slope of the gluon and quark fragmentation functions are compared in dependence of 
xe superimposed by the result of the DGLAP fits. The shaded area indicates the statistical uncertainty of Ca- 
The data points were obtained from power law fits to each xe interval individually. As expected from the structure 
of the DGLAP equation, the scaling violations are ~ 2 times larger for gluons than for quarks. Furthermore, the 
observed scaling violation in quark jets is in very good agreement with the measurements of the TASSO Collaboration 
[|| (already visible in Fig. ||). 



C. Scale dependence of the gluon and quark jet multiplicities 

There is a long standing QCD prediction that in the limit of large parton (or jet) energies the ratio r n = (N g )/(N q ) 
of the gluon over the quark multiplicity should be equal to Ca/Cf 0- Following a NNLO calculation |l|, r n is 
reduced by ~ 10% and nearly independent of energy. Fig. |8| shows the scale dependence of the jet multiplicity 
and Fig. ^| shows this ratio as a function of the hardness kh- Obviously, the ratio is much lower than expected; 
consequently neither the LO nor the NNLO prediction for r n can describe the measured data. 

To describe the quark and gluon jet multiplicities, one would set (N q ) = (N pert } as calculated from some QCD 



approach 14 la] and (N g ) = r„ ■ (N pert ), with t(kh) = Ca/Cf[1 — ?"i7o(as) — r 2~/o( a s)}- A modification of this 



ansatz had to be performed to obtain a sensible description of the measured data: 

(N q ){ KH ) = (N pert )(K H ) + K (2) 
(N g ){ KH ) = (N pert )( KH ) ■ r n + N° (3) 

The introduction of the phenomcnological offsets Nq and Nq which are assumed to be constant, accounts for expected 
differences of the fragmentation of the leading quark or gluon. 
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FIG. 8. Average charged multiplicity for light quark and 
gluon jets as a function of kh fitted with EqnJ| and ^| 
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FIG. 9. Ratio of gluon and quark jet multiplicities and 
their slopes as a function of kh 



The quality of the fit increases from x 2 /ndf—14.5 to 0.85 when introducing these additional terms. The difference 
N$-Nq yields about 2, consistent with about one additional instable primary particle built in the primary quark 
fragmentation. Unlike the multiplicity ratio itself, the ratio of derivatives of the gluon and quark multiplicities with 
scale d(N g ){KH) / Okh / d(N q ) (k^/Okh is in the predicted range of ~ 2 already at very small scale values (see Fig. ^). 
This corresponds to the expectation that if the multiplicity ratio is Ca/Cf in the limit of very large energies, the 
ratio of the slopes shows the same behaviour, simply due to de l'Hopital's rule. There should be less sensitivity of 
the slope ratio to non-perturbative effects than of the multiplicity ratio itself. Furthermore, the fit to the quark and 
gluon jet multiplicities extrapolates very well to the multiplicity ratios measured by Cleo (direct measurement of 
e + e~ — » T(ls) — > ggj decays) || and Opal (analysis of most energetic jets) 0, which is a confirmation of this 
analysis. 

Fitting Eq. H and || to the data yields Ca/Cf = 2.12 ± 0.10 sfat , in agreement with QCD. Nevertheless, large 
systematic errors are here to be expected for this value due to ambiguities in the assignment of particles to jets, the 
definition of the three jet region and the choice of the underlying jet scale. 



D. Measurement of Ca/Cf from three jet event multiplicities 



The uncertainties in the determination of the colour factor ratio quoted in Sec. [VC can be avoided by applying 
a MLLA prediction ]|| for three jet event multiplicities, N qqg , to the data: 

N qSg = N e+e - (2E*) + r n (pT) • |i . N e+e - (pj) - N \ (4) 



with: v\ = 2 {PqPg) ^ PqPg) : 2E* - J 7 ^, ; Pi : four - momenta 

The first term represents the multiplicity of an e + e~ event at a CMS energy of the invariant mass of the qq system, 
the second term half the multiplicity of a hypothetical e + e~ — * gg event with a CMS energy of twice the transverse 
momentum of the gluon. Coherence effects are included by the exact definition of the scales. Again, the term u -Nq" 
was introduced into the original MLLA prediction to absorb non-perturbative contributions. 
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Eq. H has been applied to strictly symmetric events with #2, O3 — (tt — #i/2) ± 1.5°. Therefore the only param- 
eter describing the event topology is the opening angle B\. N e + e -(y/s) has been chosen as N pert (y/s) JlJJl^] (as 



in Sec. IV C); the free parameters of these formulae are obtained from the multiplicities of e + e~ events as a function 
of the CMS energy. Fig. 8 shows the three jet event multiplicities as a function of the opening angle 9%. Superimposed 
is the fit with C'a/Cf and Nq as free parameters as a solid line, extrapolated to the region contaminated by two jet 
events (#1 < 35°), which is excluded from the fit. The dashed line represents the original prediction [Ca/Cf =2.25, 
and Nq = 0). The need for some kind of correction to account for non-perturbative effects is obvious, and the ansatz 
of adding a constant term is a simple but effective choice. The fit yields a very precise accurate measurement of 
Ca/Cf ■ 



-- = 2.246 ± 0.062 stQt ± 0.080 sas ± 0.095 t/ieo . 
Of 
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The prediction of the multiplicity ratio given in 
flifl has been tried as an alternative to the given 
r n . Although this calculation takes recoil effects 
into account, a non-perturbative offset term is 
still required. The prediction differs by about 
10% from fl|] in the NNLO term. As it does 
not reproduce the colour factor ratio contained 
in the fragmentation models which describe the 
data well, it has not been applied in this analysis. 
Fig. [To] shows the obtained result for the colour 
factor ratio in comparison of those obtained from 
four jet analyses. As far as we know our result is 
the most accurate measurement of Ca/Cf so far. 

In order to illustrate comprehensively the con- 
tents of the measurement of the three jet mul- 
tiplicity we compare in Fig. |ll| the multiplicity 
corresponding to a gg and a qq final state. The qq 
multiplicity is taken to be the multiplicity mea- 
sured in e + e~ annihilation corrected for the bb contribution. The gg multiplicity at low scale values is taken from 
the CLEO measurement ||, for which no systematic error was specified. At higher scale, twice the difference of the 
three jet multiplicity and the qq term (the first term in Eqn. ^) is interpreted as the gg multiplicity. The dashed 
curve through the qq points is a fit of N pert (y~S). The gg line is the perturbative expectation for back-to-back gluons 
according to the second term of Eqn. ||. The plot shows again that the increase of the gg multiplicity with scale is 
about twice as big as in the qq case, illustrating the large gluon-to-quark colour factor ratio Ca/Cf ■ 
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FIG. 9. Three jet event multiplicities 
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FIG. 10. The colour factor group plot 
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FIG. 11. Comparison of the charged hadron multiplicity 
for an initial qq and a gg pair as function of the scale. The 
dashed curve is a fit according of N per t{^/S), the full line is 
twice the second term of Eqn. ^. The shaded band indicates 
the uncertainty due to the error of Nq. 
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